Accuracy of direct gradient sensing by single cells 
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Many types of cells are able to accurately sense shallow gradients of chemicals across their diam- 
eters, allowing the cells to move towards or away from chemical sources. This chemotactic ability 
relies on the remarkable capacity of cells to infer gradients from particles randomly arriving at cell- 
surface receptors by diffusion. Whereas the physical limits of concentration sensing by cells have 
been explored, there is no theory for the physical limits of gradient sensing. Here, we derive such 
a theory, using as models a perfectly absorbing sphere and a perfectly monitoring sphere, which, 
respectively, infer gradients from the absorbed surface particle density or the positions of freely 
diffusing particles inside a spherical volume. We find that the perfectly absorbing sphere is superior 
to the perfectly monitoring sphere, both for concentration and gradient sensing, since previously 
observed particles are never remeasured. The superiority of the absorbing sphere helps explain 
the presence at the surfaces of cells of signal degrading enzymes, such as PDE for cAMP in Dtc- 
tyostelium discoideum (Dicty) and BARl for mating factor a in Saccharomyces cerevisiae (budding 
yeast). Quantitatively, our theory compares favorably to recent measurements of Dicty moving up 
a cAMP gradient, suggesting these cells operate near the physical limits of gradient detection. 



Cells are able to sense gradients of chemical concen- 
tration with extremely high sensitivity. This is done ei- 
ther directly, by measuring spatial gradients across the 
cell diameter, or indirectly, by temporally sensing gradi- 
ents while moving. In temporal sensing, a cell modifies 
its swimming behavior according to whether a chemical 
concentration is rising or falling in time [l| . This mode of 
sensing is typical of small, fast moving bacteria such as 
Escherichia coli, which can respond to changes in concen- 
tration as low as 3.2 nM of the attractant aspartate [3. In 
contrast, direct spatial sensing is prevalent among larger, 
single-celled eukaryotic organisms such as the slime mold 
Dictyostelium discoideum (Dicty) and the yeast Saccha- 
romyces cerevisiae Dicty cells are able to sense 
a concentration difference of only 1-5% across the cell 
0, corresponding to a difference in receptor occupancy 
between front and back of only 5 receptors @. Spatial 
sensing is also performed with high accuracy by cells of 
the immune system including neutrophils and lympho- 
cytes 01, as well as by growing synaptic cells and tumor 
cells. While there has been great progress in understand- 
ing the limits of concentration sensing and signaling in 
bacteria such as E. coli d, i, [13, IHl, 113, El, very lit- 
tle is know about the theoretical limits of direct gradient 
sensing by eukaryotic cells. 

In a recent set of experiments, van Haastert and 
Postma [6] measured the Chemotactic Index of Dicty 
cells in a cAMP gradient (Fig. 1, symbols). Chemo- 
tactic Index is defined as the distance moved in the di- 
rection of the cAMP gradient divided by the total dis- 
tance moved. To obtain the data in Fig. 1, van Haastert 
and Postma used a pipette containing different concen- 



trations of cAMP. Diffusion of cAMP out of the pipette 
established a steady-state cAMP gradient, with magni- 
tude a function of distance from the pipette. Chemotaxis 
was observed for cells as far as 700 /im from a pipette 
filled with 10~* M cAMP, corresponding to a mean con- 
centration of 7 nM and a gradient of only 0.01 nM//ini. 
This remarkable chemotactic ability raises the question - 
how closely does gradient sensing by Dicty cells compare 
with the fundamental limits on gradient sensing set by 
diffusion? 

Here we derive the fundamental limits of gradient sens- 
ing using two models for cells: a perfectly absorbing 
sphere and a perfectly monitoring sphere. Within the 
theory, gradients are estimated by comparing the discrete 
positions of particles, either absorbed on the surface of 
the sphere or measured inside the spherical volume, with 
the expected continuous distribution originating from a 
particular gradient. We find that a perfectly absorb- 
ing sphere is superior to a perfectly monitoring sphere 
for sensing both concentrations and gradients (Table I), 
since previously observed particles are never remeasured. 
Quantitatively, our theory (Fig. 1, solid curves) com- 
pares favorably with recent measurements of Dicty cells 
migrating to a cAMP-fiUed pipette suggesting that 
chemotactic ability of Dicty approaches the fundamental 
limits set by diffusion. 

I. LIMITS OF CONCENTRATION SENSING 

In this section, we consider the limits of concentration 
sensing set by particle diffusion. Consider as a measure- 
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ment device a spherical cell of radius a that can measure 
the local concentration of a certain dissolved chemical. 
Such an idealized device may make measurements 
following two different strategies: (1) The device can 
either act as a perfectly absorbing sphere and record the 
number of absorbed particles on its surface or (2) act 
as a perfectly monitoring sphere and count the number 
of particles inside its volume. In either case, from 
the number of particles, an estimate of the chemical 
concentration can be obtained. However, these estimates 
have an intrinsic uncertainty due to the randomness of 
particle diffusion. 

Perfectly Absorbing Sphere. For the perfectly ab- 
sorbing sphere the uncertainty in measuring a back- 
ground chemical concentration c is straightforward to 
derive. At steady state, the average particle current 
impinging on the sphere is J = 47r_Dac, where D is 
the chemical diffusion constant. The average number of 
particles absorbed in time T is TV = AwDacT. Since 
the particles are independent, N is Poisson distributed, 
i.e. {{6N)^) = (N). Therefore, the perfectly absorbing 
sphere has a concentration-measurement uncertainty of 
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Perfectly Monitoring Sphere. The perfectly moni- 
toring sphere was introduced by Berg and Purcell [S] as 
a parameter-free model for a cell that "perfectly" binds 
and releases all ligands that contact its surface. To quan- 
tify the time a diffusing particle spends in the cell's vicin- 
ity and is therefore capable of being measured. Berg and 
Purcell treated the cell as a permeable sphere that in- 
fers the particle concentration by counting the number of 
particles N inside its volume, and improves accuracy by 
averaging over several statistically independent measure- 
ments. A simple estimate for the resulting uncertainty in 
concentration can be obtained as follows: the number N 
is Poisson distributed and the cell counts appoximatcly 
N = a^c particles in its volume at any time. During a 
time r, the cell can make A^meas = T/{a? /D) indepen- 
dent measurements, where a^/D is the typical turnover 
time for particles inside the sphere, leading to 
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Berg and Purcell derived the exact concentration- 
measurement uncertainty for a perfectly monitoring 
sphere ( "perfect instrument" ) from the time correlations 
of particles inside the sphere, and obtained 
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which is identical to the estimate in Eq. up to a nu- 



FIG. 1: Comparison between the Chemotactic Index deter- 
mined by experiment (symbols and dashed lines) and our 
theory (solid curves). Chemotactic Index is defined as the 
distance moved by a cell in the direction of a gradient di- 
vided by the total distance moved. Experimental data was 
obtained by van Haastert and Postma Q] from DictyosteUum 
discoideum cells migrating toward pipettes containing four 
different cAMP concentrations, 0.1 (squares), 1 /xM (tri- 
angles), 10 (diamonds), 100 /xM (circles). The theoretical 
curves were obtained for a perfectly absorbing sphere using a 
single fitting parameter Da^T = 1.2 ■ 10^ Mm^i corresponding 
to, e.g., a cAMP diffusion constant of D = 300 /xm'^/s, a cell 
radius of a = 5 /xm, and an averaging time T = 3.2 s, using 
the gradient profiles from Ref. 0] and the Chemotactic Index 
from Eq. [511 Experimentally, the Chemotactic Index only 
reaches approximately 0.9 at zero distance, so we rescale our 
theory curves by 0.9. Inset: Chemotactic Index as a function 
of (c°)^/co in units of nM/(/im)^, where and Co are the 
gradient and concentration, respectively. 



merical prefactor. 

However, notice that the concentration-measurement 
uncertainty of the perfectly absorbing sphere is actually 
smaller than that of a perfectly monitoring sphere of 
the same size, because the perfectly absorbing sphere re- 
moves particles from the environment, and hence does 
not measure the same particle more than once. 



II. LIMITS OF GRADIENT SENSING 

Now consider the perfectly absorbing sphere and the 
perfectly monitoring sphere as devices for measuring 
the local gradient of a certain dissolved chemical. In 
both cases, measurements of discrete particles can be 
compared with the expected continuous distribution of 
particles originating from a particular gradient (Fig. 
2a), and hence, the gradient can be estimated. Here, we 
present in brief a theoretical derivation of the intrinsic 
uncertainty of gradient sensing (for details see support- 
ing information). We find that the intrinsic uncertainty 
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FIG. 2: Idealized models for gradient sensing by a cell. The 
gradient points along the z-axis, which is shown horizontally, 
(a) Continuum model for a perfectly absorbing sphere. The 
mean particle current density j{9) impinging on the sphere 
has axial symmetry; 9 measures the angle with respect to 
the z-axis. At steady-state, the particle concentration c is 
zero immediately outside the perfectly absorbing sphere, as 
shown schematically by the red curve superposed on the dot- 
ted background gradient, (b) Discrete particle model for the 
perfectly absorbing sphere. From the number and positions 
of particles absorbed during time T, the background particle 
concentration and gradient can be estimated, (c) Perfectly 
monitoring sphere. Particles diffuse in and out of the sphere 
without resistance. By monitoring, for a time T, the number 
and positions of particles inside the sphere, the background 
concentration and gradient can be estimated. 



is independent of the actual gradient present, and is 
always much smaller (by a factor of 7/60 ~ 12%) for the 
perfectly absorbing sphere. 



Perfectly Absorbing Sphere. The average particle 
current density j = —DVc impinging on the surface of 
a perfectly absorbing sphere of radius a at steady state 
follows from the diffusion equation, V^c = 0, and is given 
in polar coordinates by 

j{e,cP) = ^+3Dcr-e{9,cP), (4) 
a 

where cq is a constant background concentration, 
Cr is the background gradient, and e{0,(f>) = 
(cos sin 6, sin (p sin 9, cos 6) is a unit vector normal to the 
surface of the sphere (see Fig. 2a) . To best estimate the 
chemical gradient from an observed discrete density of 
particles absorbed at the surface of the sphere during 
time T (Fig. 2b), requires fitting the observed density 
= ^i)' where N is the total number of 

absorbed particles, to the expected density j{9, (j)) T from 
Eq. [31 Since the estimates of the components of the 
gradient in the x,y^ and z directions are independent, 
without loss of generality, we consider only the gradient 
estimate in the z-direction, i.e. = dc/dz, and later 
generalize to an arbitrary gradient. From the best fit, 
the estimate for the gradient in the z-direction after ab- 
sorption of particles for a time T is given by 

^ / a°^^cos9dA ^ 'EZiCos9, 

where 9i is the polar angle of the ith absorbed particle. 
We are interested in the uncertainty (accuracy) of the 
gradient measurement, which is given by the variance 

{{Sc.?) = {cD - {c.? (6) 
(E»=i cos^ Si) + (Ejli E^^j cos COS 

(47rL>a2r)2 

^ (Eticos^^.) _ (jV)(cos2g) _ CO 

{AT:Da^TY {AirDa^TY UnDa^T' 

The derivation of Eq. [S] made use of the indepen- 
dence of the particles to factorize the expectation value 
as (EtiE.'^.cos^.cos^,) = {N{N - 1)) {cos 9 f = 
(A^) 2 (cos 0)2, since the number of absorbed particles N 
is Poisson distributed. We also used (N) — AirDacoT, as 
well as {coa'^9) = 1/3. (The relation {cos'^ 9,) = (cos^ 6*) 
for absorbed particles holds even in the presence of a true 
gradient in the z direction since the gradient-weighted 
contribution {cos^ 9) is zero.) 



In three dimensions, the total uncertainty of the gra- 
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dient, normalized by co/a, is given by 



ini are therefore given by 
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with the factor of 3 arising because each component 
of the gradient contributes independently to the total 
uncertainty. This result for the uncertainty in gradient 
sensing is independent of the magnitude of the actual 
gradient present, including the case when no actual 
gradient is present. Curiously, the result is numerically 
identical to the concentration-measurement uncertainty 

(Eq. m- 



Perfectly Monitoring Sphere. 

Here, we extend Berg and Purcell's analysis of the per- 
fectly monitoring sphere ("perfect instrument") to in- 
clude gradient sensing. Specifically, we assume that the 
monitoring sphere measures not only the number but 
also the positions of all particles in its volume (Fig. 2c). 
The best estimate of the gradient is obtained by fitting 
a concentration gradient with a c = co + Cr ■ f to the 
observed time-averaged number density ^ / dtpohs {t) — 
T obtained by measuring the ex- 

act positions of all the particles inside the volume of the 
sphere for a time T. As above we focus on one component 
of the gradient, namely the gradient in the ^-direction Cz, 
and obtain as a best estimate 
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where the integral JdV is over the volume of the sphere. 
We are interested in the variance of this estimated gra- 
dient 
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where we have used JdVz^ = 47ra^/15 and have defined 
mz,T JdtJ dVzpohsit), (10) 

namely mz,T is the time-averaged total z-coordinate of 
particles inside the sphere. The expectation values in Eq. 
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- dt{m,{t)), 

where the quantity (t) is the total of the z-coordinates 
of all the particles inside the sphere at time t. To cal- 
culate mz{t), we consider the sphere embedded inside a 
much larger volume containing a total of AI particles. 
Then, niz{t) = X)f=i-^«(^)' where Zi is the z-coordinate 
of particle i if this particle is inside the sphere and is zero 
otherwise. On average there will be (N) = ^na^co parti- 
cles inside the sphere at any time t. The auto-correlation 
function {mz(t)mz{t')) of particles inside the sphere at 
time t and time t' can consequently be calculated as 

M M 

{mz{t)mz{t')) = (^^z,(i)z,(t')) (12) 

= ^TTa^cou{t-t') + {mz{t))^, 

where we have defined u{t — t') — {z(t)z{t')) for a single 
particle. Substituting Eqs. [Til and [T^ into Eq. [2] results 
in 

{(SCz) ) = 



Ana^T^Jo '^^ Jo ^^3) 



By defining a correlation time Tz ~ (l/a^) J^dTu{T) 
for the coordinate z{t), the double time integral in Eq. 
[I3]can be simplified, provided the time T is much larger 
than Tz- Using time-reversal symmetry u(t) = u{—t) 
for equilibrium diffusion (assuming small gradients), the 
variance simplifies to 



{{^Czf) = 



(14) 



The remaining task is to calculate r^, the probability 
that a particle with coordinate z inside the sphere at time 
t = is still (or again) inside the sphere at a later time 
T. We first consider the case in which the background 
chemical concentration is uniform, and later consider the 
presence of an actual gradient. Based on the solution 
of the diffusion equation in three dimensions, if a unit 
amount of chemical is released at point r', the concen- 
tration at point r at a later time r is given by 
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Using the result for the time integral from Ref. 

= 1 
dr f{\r- f"\,i^ 
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the correlation time can be expressed as a volume inte- 
gral over the sphere (the initial coordinate r ' is uniform 
in the sphere because we assume a uniform background 
chemical concentration) 
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dT dVz dV z'f{\r-f"\,T){17) 



dV z / dV' 



where r = |r]. The function 
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is analogous to the potential of a charge density in elec- 
trostatics, specifically the charge density p{z') = z' = 
p(r \ 0') = r' cos Q' for r' < a and p{r ') = for r' > a. To 
solve the final integral in Eq. 1171 we perform a multipole 
expansion of 4^(r, 9) in terms of Legendre polynomials 
P/(cos0), exploiting the rotational symmetry of ^(r,9) 
about the z axis [l3|, leading to 
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Now consider the contribution from an additional gra- 
dient with zero mean over the volume of the sphere. We 
need to perform the integrals over a non-uniform dis- 
tribution in Eq. 1171 Since a gradient along the z-axis 
contributes a factor cos 9, which leads to the vanishing 
integral sin^cos"^ 9d9 = 0, we conclude that only the 
constant background contributes to the uncertainty in 
the gradient measurement. Therefore, Eq. [19] for re- 
mains true even when an actual gradient is present. 

The result for in Eq. [11] can be used in Eq. [T3]to ob- 
tain the normalized uncertainty of gradient measurement 
by the perfectly monitoring sphere 
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where all three components of the gradient contribute 
independently. Hence, the perfectly monitoring sphere 
is not only inferior to the perfectly absorbing sphere for 
concentration sensing by a factor of 12/5 in variance 
(c/. Eqs. [l]and[3]), but is also inferior by an even larger 
factor of 60/7 for gradient sensing (c/. Eq. [7]). 



III. COMPARISON WITH EXPERIMENT 

Van Haastert and Postma Q recently measured the 
Chemotactic Index of Dicty cells in a cAMP gradient 
[f| (Fig. 1, symbols). They used a pipette containing 
different concentrations of cAMP to establish a distance- 
dependent steady-state cAMP gradients. The Chemo- 
tactic Index was defined as the distance moved by the 
cell in the direction of the cAMP gradient divided by the 
total distance moved. How does the observed chemotac- 
tic ability of Dicty compare with the fundamental limits 
on gradient sensing set by diffusion? To facilitate com- 
parison to the results of van Haastert and Postma Q , we 
have calculated the optimal Chemotactic Index for a cell 
acting as a perfectly absorbing sphere. 

To obtain the optimal Chemotactic Index CI, we as- 
sume that after averaging for a time T, a cell moves at 
a constant velocity in the direction of the estimated gra- 
dient. If we take the actual gradient to point in the z- 
direction, then the chemotactic index for one run i is 
simply cos 9i , where 9i is the angle between the true gra- 
dient and the estimated gradient. If the velocity and run 
time are the same for each run, leading to a constant run 
length /, then the average Chemotactic Index is given by 



CI = 
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To evaluate {cos 9i) for a perfectly absorbing sphere, 
we use our result (Eq. [7]) for the variance of the es- 
timated gradient in each direction, e.g. {{5cx,y,zY) — 
co/(127r Z?a'^T). Assuming a Gaussian distribution with 
these variances, as well as an actual gradient with mean 
value c° in the z direction, the 2-dimensional distribution 
of estimated gradients cp = [cx, Cz) is given by 



Pco{Cx,Cz) 



1 f,-cl/{2al)-{c,-c<if/C2'yl)^ (22) 



where ax,z = \/ ((<5cx,z)2). From this distribution, we can 
obtain the optimal Chemotactic Index 



CI - (cos0) = (cjy4T^) 
^e-y[I,{y)+h{y)l 



(23) 



where y = (c°)2/(4(t^) and /o(i) is the first (second) order 
modified Bessel function of the first kind. Fig. 1 shows a 
comparison of the optimal CI (solid curves) with the data 
of Ref. ^ . Importantly, the comparison relies on only 
a single global fitting parameter representing gradient- 
sensing ability, namely the product Do?T where D is the 
diffusion constant, a is the cell diameter, and T is the 
averaging time. Based on the estimates D = 300/zm2/s 
and a = 5iim Q, the averaging time is predicted to be 
about T = 3.2 s. (The perfect-monitor model yields an 
identical curve, but with a longer inferred averaging time 
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Ratio absorber/monitor 


Concentration: ^^^"^^ ^ 

2° 

Gradient: i'^^yK} 

(co/a)2 


4TrDacoT ^ 
1 

iirDacoT 


3 

b-rrDacoT 
15 

r-rrDacoT 


2.4 

f ^ 8.6 



TABLE I: Uncertainties in measured concentration and concentration gradient for two idealized cell models: a perfectly absorb- 
ing sphere (second column) and a perfectly monitoring sphere (third column). Also provided is the ratio of the uncertainties 
of the absorber and monitor. Parameters: diffusion constant D, radius of sphere o, averaging time T, and average chemical 
concentration cq. 



T = 27.5 s.) The theory for the optimal Chemotactic In- 
dex matches the experiment rather well. Eq. [53] further 
predicts that the Chemotactic Index depends on the gra- 
dient and the concentration cq only through the com- 
bination (c5)^/co (Fig. 1, inset). Intuitively, a/ (c°)^/co 
measures the signal to noise ratio - the signal is pro- 
portional to the true gradient |c^|, while the noise from 
particle diffusion scales as y/co- More generally, the op- 
timal Chemotactic Index depends on all variables only 
through the combination {c^)^Da^T/co. Moreover, the 
theory predicts the full distribution of run angles (see in- 
set of Fig. SI in supporting information), which can be 
obtained by integrating Eq. [22] in the radial direction for 
each angle 0. 



IV. DISCUSSION 

Many types of cells arc known to measure spatial chem- 
ical gradients directly with high accuracy. In particu- 
lar, Dictyostelium discoideum (Dicty) is well known to 
measure extremely shallow cAMP gradients important 
for fruiting body formation [1, @| and Saccharomyces 
cerevisiae (budding yeast) detects shallow gradients of 
a mating pheromone [Tsj . Accurate spatial sensing is 
also performed by cells of the immune system includ- 
ing neutrophils and lymphocytes ■ The question arises 
what are the fundamental limits of gradient sensing set 
by chemical diffusion? Here we derived such limits using 
as model cells a perfectly absorbing sphere and a per- 
fectly monitoring sphere Q . Within the theory, gradients 
are estimated by comparing the discrete distribution of 
observed locations of particles, either absorbed on the 
surface of the sphere (Fig. 2b) or measured inside the 
sphere (Fig. 2c), with the expected continuous distribu- 
tion originating from a gradient (Fig. 2a). We find that a 
perfectly absorbing sphere is superior to a perfectly mon- 
itoring sphere for concentration and gradient sensing by 
respective factors of 12/5 (= 2.4) and 60/7 (« 8.6) (Table 
I), since the perfectly absorbing sphere prevents rebind- 
ing of already measured particles. Indeed, the results 
presented here for the perfectly absorbing sphere repre- 
sent the true fundamental limits of both concentration 
and gradient sensing by cells. Our theory for the limits 
of gradient sensing compares favorably with recent mea- 
surements by van Haastert and Postma of Dicty cells 
migrating to a cAMP-fiUcd pipette (Fig. 1), suggesting 
that Dicty chemotaxis approaches the fundamental limits 



set by cAMP diffusion. 

The marked superiority of the perfect absorber for con- 
centration and gradient sensing leads us to conjecture 
that cells may have developed mechanisms to absorb lig- 
ands so as to prevent their rebinding. Such absorption 
could be implemented by ligand or ligand-receptor in- 
ternalization or by degradation of bound ligands. Even 
degradation of ligands without measurement could be ad- 
vantageous. For example, a perfect absorber that mea- 
sures only a fraction / of incident particles has the same 
uncertainties given in Table I but with an effective mea- 
surement time Toff = /T. Such an absorbing cell that 
measures only 12% of absorbed particles can still mea- 
sure gradients as accurately as a perfectly monitoring 
sphere. Similarly, an absorbing sphere of radius a with 
two small measurement patches at its poles of radius 
s (s « a), i.e. with a measuring surface-area frac- 
tion /{2a?'), effectively reduces the averaging time for 
polc-to-pole gradients to 3s^T/(2a^) (supporting infor- 
mation). Consequently, a measuring surface fraction as 
small as 4% yields the same uncertainty as the monitor- 
ing sphere. 

In fact, there are numerous examples in biology of 
ligand-receptor internalization [l^ and ligand degrada- 
tion on cell surfaces, which we speculate, might be 
related to gradient sensing. (1) Although many G- 
protein-coupled receptors are internalized by endocyto- 
sis the cAMP receptor cARl in Dicty is not [ll |. 
However, Dicty produces two forms of cyclic nucleotide 
phosphodiesterase (PDE) , which degrade external cAMP 
[1^ |20| . [2l| . One form is membrane-bound (mPDE) 
and effectively turns Dicty into an absorber, whereas the 
other form is soluble (ePDE). The membrane-bound form 
mPDE only accumulates during cell aggregration, sup- 
porting the idea that degradation of cAMP at the mem- 
brane helps accurate gradient sensing and navigation. 
Indeed, cells lacking mPDE display cell-autonomous 
chemotaxis defects even in mixed aggregates with iso- 
genic wild- type cells [2l|. Interestingly, there is good 
evidence that Dicty cells do carry out G-protein-coupled 
receptor mediated endocytosis of folic acid, another ma- 
jor Dicty chemoattractant 2lM- (2) In budding yeast, 
the receptor Ste2 binds a-factor pheromone, initiating a 
mating response including directed growth ( "shmooing" ) 
towards a potential mating partner. Ligand-bound Ste2 
undergoes internalization by endocytosis [23(. Further- 
more, the protease Barl degrades a pheromone exter- 
nally [U I25I , and may be largely membrane associated 
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[26|. (3) There are many examples of ligand-receptor 
internalization in developmental biology. For example, 
primordial germ cells in zebrafish migrate towards the 
chemokine SDF-la that activates the receptor CXCR4b. 
Ligand-induced CXCR4b internalization is required for 
precise arrival of germ cells at their target destination 
|27| . These examples suggest a correlation between lig- 
and internalization/degradation and the accuracy of cell 
polarization and movement. In presenting these admit- 
tedly speculative examples, our hope is to raise interest, 
across fields, in how the constraints of gradient-sensing 
accuracy may have shaped cellular sensing systems. 

While the absorption of ligands can improve gradient 
sensing, there is an inherent problem for an absorbing cell 
to measure a gradient while moving. An absorbing cell 
moving in a uniform concentration creates an apparent 
gradient due to an increased flux of incoming particles 
at its front and a decreased flux of particles at its back 
Q . Using the model of a spherical cell, the ratio of fluxes 
between front and back hemispheres is given by 



R=l 



3avQ 

IT' 



(24) 



where wq is the cell velocity, a is the cell radius, and D is 
the particle diffusion constant. On the other hand, the 
flux ratio of a stationary spherical cell in a gradient |Vc| 
with uniform background concentration background Co is 
given by 



R = 1 



3a|Vc| 
Co 



Hence, a moving cell sees an apparent gradient 
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As an example, chemotaxis of Dicty to cAMP is observed 
at a mean concentration of 7nM in a gradient of only 0.01 
nM//im [6]. A Dicty cell moving with a typical speed of 
0.2 /im/s at the same mean concentration but without a 
gradient creates an apparent gradient of about half the 
real gradient. There are several ways out of this dilemma. 
(1) Cells could separate measurement from movement at 
low gradients, e.g. by stopping, measuring the gradient, 
and then moving. Dicty would only need to stop for 
about a? /D « 0.1s based on cell radius of a = 5/zm and 
diffusion constant D = BOO/zm^/s. (2) Cells could sense 
gradients transverse to their direction of motion. This is 
particularly advantageous for fast moving cells {e.g. bac- 
teria) for which the apparent gradient can become more 
than 100 times steeper than the actual gradient @. (In- 
deed, the oxygen-sensing marine bacterium Thiovulum 
majus directly senses gradients transverse to its direc- 
tion of motion [1^.) Interestingly, Dicty cells moving on 
agar in the absence of a gradient appear to combine these 
two strategies. Qualitatively, the tips of elongated mov- 
ing cells slow down and ffatten, often producing two or 



more distinct pseudopods. Cells then elongate and move 
(~ one cell length) in the direction of one of the pseu- 
dopods before the process is repeated (Liang Li and Ted 
Cox, personal communication). By this strategy, Dicty 
may avoid locking onto a false, movement-generated gra- 
dient. (3) In principle, cells could compensate for the 
apparent motion-generated gradient, cither by internal 
signal processing or by external chemical secretion. In 
fact, Dicty cell do secrete cAMR primarily from their 
trailing edge during movement [Jl, but this cAMP se- 
cretion serves dominantly to facilitate cell aggregation, 
including cells following cells during streaming. Given 
the complex role of cAMP in Dicty aggregation, studies 
of Dicty chemotaxis using gradients of folate, which is 
absorbed [12] but apparently not secreted by Dicty, may 
ultimately prove simpler to interpret. 

Our models of the absorbing and the monitoring 
spheres neglect all biochemical reactions, such as particle- 
receptor binding and downstream signaling, which could 
significantly increase measurement uncertainty beyond 
the fundamental limits described here. To study the ef- 
fects of particle-receptor binding, we extended a formal- 
ism for the uncertainty of concentration sensing, recently 
developed by Bialek and Setayeshgar [l3|, to gradient 
sensing. We found that the measurement uncertainty al- 
lowing ligand rebinding is larger than the measurement 
uncertainty without rebinding, confirming the superiority 
of the absorber over the monitor (details will be published 
elsewhere). A number of mechanistic models for gradi- 
ent sensing and chemotaxis have addressed the impor- 
tant questions of c ell p olarization, signal amplification, 
and adaptation S^, 3ll [3^ . [H, [s^, [H, [sl] , cell movement 
of individual cells |3 7| . [S q. cell aggregation with cAMP 
degradation by PD E |39| . as well as sensing of fluctuating 
concentrations [l^, l40ll4l| . Our results on the fundamen- 
tal limits of gradient sensing complement these models, 
and may ultimately help lead to a comprehensive descrip- 
tion of eukaryotic chemotaxis [i^l ■ 

Finally, we remark that the experiments by van 
Haastert and Postma used stationary spatial gradients 
0. Cells in such gradients might profit from remember- 
ing their direction of motion 43] , and evidence for such 
internal memory was recently obtained [H, [U, S^- It 
therefore might prove interesting to measure the chemo- 
tactic index for randomly changing gradients, to find out 
if cells indeed use their memory to improve chemotaxis. 
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I. ACCURACY OF CONCENTRATION AND GRADIENT SENSING BY A PER- 
FECTLY ABSORBING SPHERE 



Consider as a measurement device a sphere of radius a that is a perfect absorber for a 
certain dissolved chemical. The device can be used to measure both the local concentration 
and the local gradient of the chemical by sensing the chemical current density j impinging 
on the surface of the sphere. At steady-state, for the case of a uniform particle concentration 
Co far away from the absorbing sphere, the average current is J = AnaDco, where D is the 
particle diffusion constant Since the average number of particles N = AnDacoT absorbed 
in time T is Poisson distributed, i.e. {{SN)"^) = {N), the perfectly absorbing sphere has a 
concentration-measurement uncertainty of 

m') mr) i .... 

c2 (iV)2 AnDacT' ^ ' 

The perfectly absorbing sphere can also be used to measure the local gradient of the 
chemical. To calculate the current density j in this case, we utilize a standard analogy to 
electrostatics. In electrostatics the potential (f) and electric field ii^ in a charge-free environ- 
ment are determined by Laplace's equation, = 0, and E = — V0, respectively. As a 
result, the surface-charge density (Xchargc on a conducting sphere (boundary condition = 
at r = a) placed in an electric field of magnitude in the z-direction with an additional 
constant potential far away from the sphere is given by the superposition [2I 



1 

""charge 



-— {- -3E,cos9] (12) 



An dr 477 \^ a 

in Gaussian units, where 6 is the polar angle measured with respect to the z axis. In the 
case of a chemical-absorbing sphere, the chemical concentration c and the current density 
j obey equations analogous to the ones governing the potential and electric field E in 
electrostatics. Specifically, the spatial dependence of the concentration c follows from the 
diffusion equation at steady-state, V^c = 0, while the current density is given by j = —DVc. 
Exploiting the result from electrostatics, Eq. [T2|, the average current density impinging on 
the perfectly absorbing sphere (boundary condition c = at r = a) in a background gradient 
Cz = dc/dz in the z-direction is given by 

j(e) = ^ + 3DczCos9. (13) 



This expression can be generalized to a gradient Vc in an arbitrary direction r via 

j(^,0) = ^ + 3DVc-e(e,0), (14) 
a 

where e{6, (p) = (cos (p sin 6, sin (p sin 6, cos 6). 

To best estimate the chemical gradient from an observed density of particles absorbed at 



the surface of the sphere during time T, we fit the observed density = J2iLi ^(^ 
where N is the total number of absorbed particles, to the expected density j{6,(j))T from 
Eq. [m The best fit is obtained by minimizing the error between the observed density and 
the expected density 

Error = f [a?,"^ -C- F,™ (^^, (P)]'dA, (15) 

m=-l,0,l 

where G and Gm are the parameters to be determined. Note that in Eq. [15] we write 
the expected contribution from the gradient in terms of the spherical harmonics Yj^^{6,(j)): 
= y^sin^e"**^, Y^ = ^J^cosO, Y^ = — y^^sin^e*'^. Minimizing the error as a 
function of the parameters G and Gm (^ = —1, 0, 1) is achieved by setting d'Eiioi / dG = 
and dEiToi/ dGm = 0, which results in the best-fit values 

r a^^' dA f a^^' dA , , 



J <-yf ^ (g, (P)dA sin ee-''t> dA 

~ j\Y^-\e,<P)\^dA ~ a2 ^^^^ 



/ <^yo(^, (t))dA COS e dA 

^0 - nvo^n AM2^A ^2 (IS) 



f\Y,%6,<PWdA 



_ /<%^(g,0)rfA _ -Vl^/<^singe-'^c;A 
' ~ J\Y,\e,<f))\^dA ~ a2 ^ ^ 

and therefore the following best estimates for the background concentration and the indi- 
vidual gradient components: 

ci'G 

Co = 7^ (20) 



DT 
1 / 3 



QDT V 27r 




(G_i - Gi) (21) 



c.^^yi(G. + G-0 (22) 
^ ^Go. (23) 



6Dr V vr 



Since c^^y^z are estimates of the independent, orthogonal components of the gradient, without 
loss of generality, we consider only the gradient estimate in the 2;-direction and generalize 
to an arbitrary gradient later. From Eq. [23l the best estimate for the gradient in the 
^-direction after absorption of particles for a time T is given by 

^ /<-cosgdA ^ E^cosg. .2^. 

We are ultimately interested in the uncertainty (accuracy) of the gradient measurement 
by the absorbing sphere, which is given by the variance 

{{^c,f) = {cD - (25) 

(26) 

(27) 



(Er=i cos^ e,) + (Eti ES., COS 6, COS 6,) (Et, COS 6,^ 



(Ell cos' 0,) _ {N){cos'e) _ Co 



{A'kDo?TY {A-KDa^Tf UnDa^T' 

Deriving Eq. [27] made use of the independence of the particles to factorize the expectation 
value (EiIiE!iiCos^,cos^j) = {N{N - l))(cos^)2 = {N)'{cose)', since is Poisson 
distributed. The last step is exact. We also used (A^) = ArcDacoT, as well as (cos^6') = 1/3. 
The latter relation applied to absorbed particles holds even in the presence of a gradient in 
the z direction since the average (cos'^ 6) is zero. 

Since the gradient may point in an arbitrary direction, the total uncertainty of the gra- 
dient, normalizedby Co/a, is given by 

(co/a)2 (co/a)2 AnDacoT' ^ ' 

with the factor of 3 arising because each component of the gradient contributes inde- 
pendently to the total uncertainty. This result for the uncertainty in gradient sensing is 
independent of the magnitude of the actual gradient present (including the case when no 
gradient is present). Interestingly, the result is identical to the uncertainty in measureing 
the concentration by the same perfectly absorbing sphere (c/. Eq. [TTI) . 



Special case: absorbing sphere with two polar measurement patches 

We consider a special case in which a whole sphere of radius a is absorbing, but only two 
small circular patches of radius s (^ a) opposite to each other on the sphere measure the 



number of absorbed particles. We start with the analogues of Eqs. [TS] and [21] for the 
estimated gradient in the 2;-direction 

Go = ^ (29) 



Go _ / of' cos edA _N^-N^ ^^^^ 



2V^DT GnDs^T GnDs^T' 
where we have used cos 6' ~ 1 for patch 1 and cos^ ^ — 1 for patch 2. The variance is given 
by 

{{Sc.r) = {cD - (c.)2 (31) 



{67cDs^Ty 
Co 



(32) 
(33) 

ISnDs'^aT ^^"^^ 
where we have used (SNf) = {Nf) — (A^i)^ = {Ni) {i = 1,2) for a Poisson process, as well 
as (A^i) = tts'^DT^cq + {cz) a) /a and {N2) = ns'^DT^co — {cz) a) /a. From this we obtain the 
normalized uncertainty of the gradient in the z direction 

{{Sc.?) _ 1 2a^ 

(co/a)2 UnDacoT Ss^ ' ^ ^ 

i.e. 2a^/(3s^) times the uncertainty in the gradient using the whole sphere. Hence, the 
smaller the patch size s, the larger in the uncertainty of the gradient measurement. 

II. ACCURACY OF GRADIENT SENSING BY A PERFECTLY MONITORING 
SPHERE 

Berg and Purcell considered a "perfect instrument" for measuring the concentration of 
a certain dissolved chemical, namely a virtual sphere of radius a that could exactly count 
the number of diffusing chemical molecules inside its volume [l|. They showed that such 
an instrument, making measurements for a time T, could estimate a concentration with an 
uncertainty 

c2 " hixDacT' ^^^^ 
We note that the "perfect instrument" is actually inferior by a factor 12/5 in variance to the 

perfectly absorbing sphere for concentration measurement (Eq. [TT]) because the absorbing 



sphere removes particles from the environment, and hence does not measure the same particle 
more than once. 

We extend Berg and Purcell's analysis of the perfect instrument to include gradient 
sensing by assuming that the virtual sphere also measures the positions of all particles in 
its volume. As in the previous section, we derive a best estimate for the gradient and the 
variance of this best estimate to find the uncertainty of the gradient determination. We 
start by fitting a gradient model c — Co + f- Cr to the observed time-averaged number density 
^ J dtpohs{t) — ^ J d't^^i5{f — f{t)), obtained by measuring the exact positions of the 
particles inside the volume of the sphere for a time T. Specifically, we minimize the error 

Errors J J Pohs{t)dt - - r ■ dV. (37) 

As before we focus on one component of the gradient, namely the gradient in the 2;-direction 
Cz. By setting 9Error/9c2 = 0, we obtain as a best estimate for the z gradient 

^ J dt JdVzpohs{t) 



/ dVz^ 

We are interested in the variance of this estimated gradient 



(38) 



{{5c,f) = {cD - {c^f (39) 

/ 15 \' . ,,,,,, . X > X /,,,,, . X , V . ^ 

(40) 



V 47ra5 I ' 




dt J dVzpobsit)^ ) ~ ((/'^^/ '^^'^'^°'^^(^^) 



where we have used / dVz^ — Ana^/lb and we have defined 

dt J dVzpob.{t)^^ ^fJo 

The quantity mz{t) is the total of the z coordinates of all the particles inside the sphere 
at time t. To calculate mz{t), consider the sphere imbedded inside a much larger volume 
containing a total of M particles. This results in mz{t) = X^^i^i(t), where Zi is the z- 
coordinate of particle i if this particle is inside the sphere and zero if it is outside. If 
N{t) is the number of particles inside the sphere at time t, then on average there will be 




(A^) = Ivra^co particles inside the sphere. The auto-correlation function {mz{t)mz{t')) of 
particles inside the sphere at time t and time t' can consequently be calculated 

M M 

K(t)m.(t')) = (44) 
i=i i=i 

M MM 

= iJ^^i^Mt')) + ■ (45) 

i=l 1=1 j^i 

= {N) {z{t)z{t')) + {Nf {z(t)f (46) 
= ^vra^co u{t - t') + (m,(^))^ (47) 

where we defined u{t — t') = {z(t)z(t')). For Eq. |46] we used that ~ 
(Sj=i which is exact in the thermodynamic limit of the large embedding volume 
and large M. Furthermore, since N{t) and z{t) are independent, the factorization of the 
expectation values is exact as well. Plugging Eqs. HT] and US] into Eq. Hi] results in 



((5c.)2) = -^ / dt I dt'uoit-t'), (48) 



where the second term in Eq. HT] cancels Eq. HS] 

By introducing a correlation time = l/a^ ■ Jq°° dru^r) for z(t), the double time integral 
in Eq. HH] can be simplified, provided time T is much larger than r^. Using time-reversal 
symmetry u{t) = u{—t) for equilibrium diffusion (small gradients), the variance transforms 
into 

((^c.,^> ^ (49) 

The remaining task is to calculate r^, the probability that a particle with coordinate z inside 
the sphere at time t = is still (or again) inside the sphere at a later time r. Since the 
chemical concentration consists of a constant background concentration plus a gradient, we 
first consider the constant background and later consider the additional gradient. 

Based on the solution of the diffusion equation, if a unit amount of chemical is released 
at point r', the concentration at point r at a later time r is given by 



{AttDt 

Using the result for the time integral in Ref. [l| 



/(k"*-r1,r) = ^^— (50) 



the correlation time can be expressed as a volume integral over the sphere (the initial 
coordinate r' is uniform in the sphere because we assume a uniform background chemical 
concentration) 

= r dr [ dVz [ dV'z'f{\f-f'\,T) (52) 

3^Cf Jo J J 

= I'^'I '^'W^l = / f^^' «)■ 

where r = |f^. The function il){r,9) is analogous to the potential of a charge density that 
is rotationally symmetric around the 2;-axis in electrostatics, i.e. p{z') = z' = p{r',6') = 
r'cosO' for r' < a and p(r') = for r' > a. To solve the integral in Eq. [53l we perform a 
multipole expansion of the potential 

in terms of Legendre polynomials Pi{cos6), exploiting the rotational symmetry about the z 
axis 121. One needs to differentiate two cases: 



Case I (r' < r < a): 



1=0 ^ ^ 



^>(^>^)=>^ iSrM'Kcos^^), (55) 



where the exterior multipole moments are given 6') with p(r ') 

r' cosO'. Performing the integral, only the dipole moment survives 



yielding 

Case II (r < r' < a): 



Qi = (56) 



^>(r,^) = ^cose. (57) 



^<(r,^) = ^ VP^cos^), (58) 



1=0 



where the interior multipole moments are defined &s Ii = J dV'^jT^Pi{cos9') and, again, 
only the dipole moment survives 



yielding 

^^(r,e) = cosg. (60) 

o 

Using these expressions for the potential, the remaining integral in Eq. [53] can be performed 
by summing up the two contributions to the potential ipi^, ^) = V^<('^) ^) + i^>{^: ^) 

/ drr^ dd ?,ind co&d U^{r,e) + i)y{r,d)] = . (61) 

Jo Jo lOSL* 



Now consider the contribution from an additional gradient with zero mean over the volume 
of the sphere. We need to integrate j d^r... over a non-uniform distribution in Eqs. [53] and 
[nU Since a gradient along the 2;- axis contributes a factor cos^, leading to the vanishing 
integral ^ sin cos^ 9 d9 = 0, only the constant background contributes to the uncertainty 
in the gradient measurement. 

The result for (Eq. [6T]) can be used in Eq. [49] to obtain the normalized uncertainty of 
the gradient measurement by the perfect instrument 

((^ = 5 .62) 
(co/a)2 l-KDac^T' ^ ' 

Since each component of the gradient contributes independently, the total normalized un- 
certainty is finally 

(co/a)2 YnDacoT 

Hence, the "perfect instrument" is not only inferior to the perfectly absorbing sphere for 
concentration sensing by a factor of 12/5 in variance (c/. Eqs. [36l and [TTj) . but is also inferior 
by an even larger factor of 60/7 for gradient sensing (Eq. [28|l . 
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FIG. 1: Main panel is reproduced from main text Fig. 1, showing the Chemotactic Index of 
Dictyostelium discoideum cells migrating toward a pipette. Symbols correspond to four different 
cAMP concentrations, 0.1 fiM (squares), 1 fiM (triangles), 10 /jM (diamonds), 100 fiM (circles). 
Inset: predicted distributions of cell-movement directions at the two points of the 100 fiM cAMP 
curve indicated by arrows, distances 200 jum (solid curve and arrow) and 500 fim (dashed curve 
and arrow) from the pipette. 



